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ABSTRACT 

We perform a detailed investigation of moderate-to-high quality X-ray spectra of ten of the most luminous 
active galactic nuclei (AGNs) known at z > 4 (up to z ~ 6.28). This study includes five new XMM-Newton 
observations and five archived X-ray observations (four by XMM-Newton and one by Chandra). We find that 
the X-ray power-law photon indices of our sample, composed of eight radio-quiet sources and two that are 
moderately radio loud, are not significantly different from those of lower redshift AGNs. The upper limits 



obtained on intrinsic neutral hydrogen column densities, A/h 10 -10 cm , indicate that these AGNs are 
not significantly absorbed. A joint fit performed on our eight radio-quiet sources, with a total of ^7000 pho- 
tons, constrains the mean photon index of z > 4 radio-quiet AGNs to T=1.97 
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with no detectable intrinsic 



dispersion from source to source. We also obtain a strong constraint on the mean intrinsic column density, 



A/h 3x 10 cm , showing that optically selected radio-quiet AGNs at z > 4 are, on average, not more ab- 
sorbed than their lower-redshift counterparts. All this suggests that the X-ray production mechanism and the 
central environment in radio-quiet AGNs have not significantly evolved over cosmic time. The mean equivalent 
width of a putative neutral narrow Fe Ka line is constrained to be ,$,190 eV, and similarly we place constraints 
on the mean Compton reflection component (R < 1.2). None of the AGNs varied on short (~1 hr) timescales, 
but on longer timescales (months-to-years) strong variability is observed in four of the sources. In particular, 
the X-ray flux of the z=5.41 radio-quiet AGN SDSS 0231-0728 dropped by a factor of ~4 over a rest-frame 
period of 73 d. This is the most extreme X-ray variation observed in a luminous z > 4 radio-quiet AGN. 
Subject headings: galaxies: active - galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

Active galactic nuclei (AGNs) at z > 4 are valuable cosmo- 
logical probes of the physical environment in the ~1 Gyr old 
Universe, following the re-ionization epoch at the end of the 
'dark ages' (e.g., Fan et al. 2002). Investigations of physi- 
cal processes in these sources, and in particular their energy 
production mechanism, are important for understanding the 
ionizing flux contribution from AGNs in the early Universe, 
black hole (BH) growth, and the evolution of mass accretion 
in galactic nuclei (see Brandt & Hasinger 2005 for a recent 
review). It is not yet clear whether the small-scale physics of 
active nuclei are sensitive to the strong large-scale environ- 
mental changes the Universe has experienced since z ~ 6. 

The penetrating X-ray emission from AGNs is of particu- 
lar interest, since it allows one to probe the innermost region 
of the accretion flow, near the central supermassive BH. Un- 
til fairly recently, X-ray detections of z > 4 AGNs were lim- 
ited to a handful of sources with poor photon statistics, pre- 
venting a detailed analysis of their X-ray spectral properties 
(Kaspi, Brandt, & Schneider 2000); this was also, in part, a 
consequence of the relatively small number ( 100) of such 
sources discovered by the year 2000. The past five years have 
seen a rapid increase in optical discoveries of z > 4 AGNs, 

1 Department of Astronomy & Astrophysics, The Pennsylvania State Uni- 
versity, University Park, PA 16802, USA; ohad@astro.psu.edu 

2 INAF - Osservatorio Astronomico di Bologna, Via Ranzani 1, 40127 
Bologna, Italy 

3 Dipartimento di Astronomia, Universita' degli Studi di Bologna, Via 
Ranzani 1, 40127 Bologna, Italy 

4 Steward Observatory, University of Arizona, 933 North Cherry Avenue, 
Tucson, AZ 85721, USA 

5 Princeton University Observatory, Peyton Hall, Princeton, NJ 08544, 
USA 



largely as a result of two surveys, the Palomar Digital Sky 
Survey (DPOSS; Djorgovski et al. 1998) and the Sloan Dig- 
ital Sky Survey (SDSS; York et al. 2000), which have raised 
the number of such sources to ^600 (e.g., Schneider et al. 
2005). A rapid increase in X-ray detections of z > 4 AGNs 
soon followed, with ^90 such sources to date. 6 Chandra and 
XMM-Newton observations of these sources have provided 
constraints on their mean global X-ray spectral properties, in 
particular their photon index and intrinsic absorption column 
density (e.g., Brandt et al. 2004; Vignali et al. 2005). How- 
ever, many of the X-ray observations are shallow, aimed at 
just detecting the AGN, and do not allow useful measurements 
of the spectral properties in individual sources (e.g., Brandt et 
al. 2002; Vignali et al. 2003a, 2005). Exposure times of sev- 
eral tens of ks with Chandra and XMM-Newton are needed to 
measure the spectral properties accurately. Such observations 
are important for investigating the energy production mecha- 
nism at this early epoch, and testing whether the X-ray spec- 
tral properties have evolved over cosmic time. 

Recent studies of X-ray spectral properties and their de- 
pendence on redshift and luminosity have led to conflicting 
conclusions. For example, while Bechtold et al. (2003) re- 
ported that the X-ray photon indices (r) of z > 4 AGNs are 
flatter (1.5±0.5) than those of nearby AGNs, Grupe et al. 
(2005) reported that their T are rather steep (2.23±0.48). 
The first study was based upon a sample of 16 radio-quiet 
sources observed by Chandra with low total counts (see Ta- 
ble 1 of Bechtold et al. 2003), and one of their objects 
(SDSS J103027.10+052455.0, hereafter SDSS 1030+0524) 
is common with the present work. The second study was 

6 See http://www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat 
for a regularly updated compilation of X-ray detections at z > 4. 
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TABLE 1 

LOG OF XMM-Newton OBSERVATIONS 



AGN 



a 

(J2000.0) 



S Galactic /V H a Obs. 

(J2000.0) (10 20 cm- 2 ) Date 



Exp. Time (ks) / Total Counts 
MOS1 MOS2 pn 



Ref. 



New Observations 



PSS0121+0347 b 


4.13 


01 21 26.2 


+03 47 06.3 


3.23 


2004 Jan 09 


25.7/139 


25.9/185 


21.7/540 




SDSS0210-0018 b 


4.77 


02 10 43.2 


-00 18 18.3 


2.66 


2005 Feb 05 


18.7/112 


18.7/117 


16.5/594 




SDSS 0231-0728 


5.41 


02 31 37.7 


-07 28 54.5 


2.90 


2004 Jan 07 


31.6/70 


31.6/101 


27.7 / 309 




PSS 0926+3055 


4.19 


09 26 36.3 


+30 55 04.9 


1.89 


2004 Nov 13 


27.8 / 289 


27.8 / 295 


23.8/1156 




PSS 1326+0743 


4.17 


13 26 11.8 


+07 43 57.5 


2.01 


2003 Dec 28 


31.3/247 


31.3/229 


27.1/963 




Archival Observations 


Q 0000-263 


4.10 


00 03 22.9 


-26 03 17 


1.67 


2002 Jun 25 


34.0 / 295 


34.0 / 322 


34.0/ 1229 


2 


BR 0351-1034 


4.35 


03 53 46.9 


-10 25 19 


4.08 


2002 Aug 23 


17.9/31 


17.9/15 


15.1 / 166 


3 


SDSS 1030+0524 


6.28 


10 30 27.1 


+05 24 55 


3.17 


2003 May 22 


44.3 / 80 


43.9/81 


43.4/333 


4 


BR 2237-0607 


4.56 


22 39 53.6 


-05 52 19 


3.84 


2003 May 17 


17.2/60 


17.3/76 


16.3 / 306 


5 



REFERENCES. — 1. This work; 2. Ferrero & Brinkman (2003); 3. Grupe et al. (2004); 4. Farrah et al. (2004); 5. Grape et al. (2005) 

"Galactic absorption column densities from Dickey & Lockman (1990), obtained using the optical coordinates of the sources given in this Table 
and using the HEASARC Nh tool at http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
'Moderately radio-loud source. 



based upon a sample of 16 sources (11 of which are ra- 
dio quiet) observed with XMM-Newton with low-to-high to- 
tal counts (see Table 1 of Grupe et al. 2005); four of those 
sources (Q 0000-263, BR 0351-1034, SDSS 1030+0524, 
and BR 2237—0607) are common with the present work, and 
we use the same X-ray data that were discussed in that study. 
The recent analysis by Vignali et al. (2005), which is based 
upon the joint X-ray spectral fitting of a sample of 48 z > 4 
radio-quiet AGNs observed with Chandra in the redshift range 
3.99-6.28, argues that the X-ray spectral properties of AGNs 
do not show any significant evolution or luminosity depen- 
dence. Our aim in this paper is to address this puzzle by uni- 
formly analyzing a sample of moderate-to-high quality X-ray 
spectra of z > 4 AGNs, with «200-1000 counts per source. 
In this paper, we consider radio-quiet and radio-loud AGNs 
separately, as radio-loud sources might exhibit X-ray emis- 
sion from their jets. The X-ray properties of the two AGN 
sub-classes differ significantly from one another; for example, 
the X-ray spectra of radio-loud AGNs are flatter than those of 
radio-quiet sources (e.g., Wilkes & Elvis 1987; Cappi et al. 
1997; Porquet et al. 2004) and, at high redshift, radio-loud 
AGNs show signs of higher intrinsic absorption (e.g., Wors- 
ley et al. 2004a, b) and larger X-ray variability (e.g., Fabian et 
al. 1999). Since radio-quiet sources represent the majority of 
luminous AGNs at all redshifts (e.g., Stern et al. 2000; Ivezic 
et al. 2002; Petric et al. 2003), our analysis focuses on the 
properties of this sub-class. 

Moderate-to-high quality spectra of only five z > 4 radio- 
quiet AGNs have been obtained with XMM-Newton and 
Chandra over the past three years, with only one such source 
having more than 1000 counts (Q 0000—263, Ferrero & 
Brinkmann 2003; see also Farrah et al. 2004; Grupe et al. 
2004; Schwartz & Virani 2004; Grupe et al. 2005). In this 
paper we almost double that number with new data from 
XMM-Newton, adding two more sources with over 1000 
counts, and study the broad-band X-ray spectral properties 
of ten AGNs with an average redshift of 4.8. In §[2] we de- 
scribe our five XMM-Newton observations, the data reduc- 
tion, and the analysis. In §[5] we expand our original sample 



by adding archived high-quality XMM-Newton and Chandra 
spectra of five z > 4 radio-quiet AGNs, and in §[4] we inves- 
tigate the mean X-ray spectral properties of our radio-quiet 
sample. Our results are discussed in § [5] and summarized in 
§ [6] Throughout the paper we use the standard cosmologi- 
cal model, with parameters Q.^=0.7, Q.m=0.3, and Hq=1Q km 
s~ l Mpc~ l (Spergel et al. 2003). 

2. NEW HIGH-QUALITY SPECTRA OF Z > 4 AGNS 
2.1. Observations and Data Reduction 

We have obtained imaging spectroscopic observations of 
five z > 4 AGNs with XMM-Newton (Jansen et al. 2001). 
Three of the sources were discovered in the DPOSS survey 
and were detected in X-rays by Vignali et al. (2003a); two 
were discovered in the SDSS survey and detected in X-rays 
by Vignali et al. (2001, 2003b). Three of the sources are ra- 
dio quiet with an upper limit on the radio flux, while two, 
PSS 0121+0347 and SDSS 0210-0018, are moderately ra- 
dio loud, following the classification of Kellermann et al. 
(1989), and have fl=300 and fl=80, respectively (Vignali et 
al. 2003a, 2001). These five sources were among the X-ray 
brightest of a sample of z > 4 AGNs previously targeted with 
Chandra snapshot ( ^ 10 ks) observations (Vignali et al. 2001, 
2003a, b). One of these sources, SDSS 0231-0728 at z=5.41 
(Anderson et al. 2001), is the highest redshift AGN detected 
with the automated SDSS detection algorithm. The short, 
higher angular resolution Chandra observations allow predic- 
tion of the number of counts expected from the XMM-Newton 
observations, and they eliminate the possibility of signifi- 
cant source confusion due to the lower angular resolution of 
XMM-Newton. No nearby contaminating sources were de- 
tected in the Chandra images of our five AGNs. 

A log of the new XMM-Newton observations appears 
in Table ^ (first five sources). The data were processed 
using standard SAS 7 V6.1.0 and FTOOLS tasks. The 
event files of all the observations were filtered to include 



7 XMM-Newton Science Analysis System. See http://xmm.vilspa.esa.es/ 
external/xmm_sw_cal/sas_frame.shtml 
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TABLE 2 

Best-Fit X-ray Spectral Parameters 



PSS 0121+0347 
8=4.13 



AGN 




r 


f E (l keV) b 


Z 2 (DOF) 


PSS 0121+0347 


<2.91 


1 S1+ 016 
1 - 81 -0.16 


7 64 +0% 
'•"^-1.01 


54.8(52) 


SDSS 0210-0018 


<4.17 


1 81+0-15 


6 45+ 78 


64.2(75) 


SDSS 0231-0728 


< 19.90 




2 10+ 059 


25.3(43) 


PSS 0926+3055 


< 1.02 


1 gg+0.08 


n.5&+\f s 


76.1(78) 


PSS 1326+0743 


<0.47 


, o 7 +0.10 


12 - 54 -?:i3 


51.2(65) 



Q 0000-263 


<0.40 


2.02; 


-0.07 
-0.07 


12.59 


+0.75 
-0.75 


71.3(82) 


BR 0351-1034 


<5.29 


1.76; 


-0.31 
-0.26 


5.24^ 


-0.96 
-0.97 


8.7(10) 


SDSS 1030+0524 


< 7.77 


1.92; 


-0.23 
-0.21 


1.39^ 


-0.23 
-0.23 


40.9(44) 


BR 2237-0607 


< 3.30 


2.09; 


-0.22 
-0.21 


3.82^ 


-0.59 
-0.57 


20.4(25) 



NOTE. — The best-fit photon index, normalization, and x 2 were obtained 
from a model consisting of a Galactic-absorbed power-law without any in- 
trinsic absorption. 

"Intrinsic column density in units of 10 22 cm -2 . Upper limits were calcu- 
lated using the intrinsically-absorbed power-law model with Galactic absorp- 
tion, and represent 90% confidence limits for each value, taking one parame- 
ter of interest (A^ 2 =2.71; e.g., Avni 1976). 

b Power-law normalization for the pn data, taken from joint fitting of all 
three EPIC detectors with the Galactic-absorbed power-law model; units are 
10~ 6 keVcm~ 2 s" 1 keV _1 . 



events with FLAG=0, and PATTERN<12 (PATTERN<4) and 
200<PI< 12000 (150<PI< 15000) for the MOS (pn) detec- 
tors. The event files of the observations of SDSS 0210-0018, 
SDSS 0231-0728, and PSS 0926+3055 were also filtered 
to remove periods of flaring activity, which were apparent in 
the light curve of the entire full-frame window. Such flar- 
ing activity, which originates in the Earth's magnetosphere, 
depends on several unpredictable factors, and it can signifi- 
cantly reduce the good-time intervals of an observation. 8 The 
observations of SDSS 0210-0018, SDSS 0231-0728, and 
PSS 0926+3055 were shortened by 16 ks, 10 ks, and 0.6 ks, 
respectively, as a result of this filtering. The event files of 
the PSS 0121+0347 and PSS 1326+0743 observations were 
not filtered in time, since the entire observation of the first 
source was performed during a moderate background flare pe- 
riod with count rates ^3 times higher than nominal values, 
and the full-frame count-rate light curve of the second source 
did not show any flaring. The exposure times listed in Table[l] 
reflect the filtered data used in the analysis. 

We searched for X-ray emission from an additional radio- 
quiet AGN, SDSS J021 102.72-000910.3 at z=4.90 (Fan et 
al. 1999), which lies close to the edge of the field of view 
in the XMM-Newton observation of SDSS 0210-0018. The 
former source is not detected, and the 3 a detection thresh- 

8 On average, ~35% of the data from XMM-Newton observations with 
exposure times greater than 10 ks are lost due to flaring (K. Kuntz, 2005, 
private communication). For more details see the XMM-Newton User 
Handbook at http://xmm.vilspa.esa.es/external/xmm_user_support/ documen- 
tation/uhb/index.html 




SDSS 0231-0728 
_=5.41 





Observed— Frame Energy (keV) 

FIG. 1. — Data, best-fit spectra, and residuals for our new XMM-Newton 
observations of five z > 4 AGNs. Open circles, filled squares, and open 
squares represent the pn, MOS1, and MOS2 data, respectively. Solid lines 
represent the best-fit model for each spectrum, and the thick line marks the 
best-fit model for the pn data. The % residuals are in units of o~ with error 
bars of size one. The inset in each panel shows 68, 90, and 99% confidence 
contours for the intrinsic absorption (N\\) and photon index (T). 

old is ^120 counts. Using PIMMS, 9 this translates to an up- 
per limit on the unabsorbed observed-frame 0.5-2 keV flux 
of ;$ 10~ 14 ergs cm~ 2 s , which is consistent with the 
3.II14X 10~ 15 ergs cm~ 2 s _1 flux measured by Vignali et al. 
(2001). 



2.2. Spectral Analysis 

To extract X-ray spectra we used a source-extraction aper- 
ture radius of 30" in the images of all three European Pho- 
ton Imaging Camera (EPIC) detectors for all but one AGN, 
PSS 0121+0347. For this source we used an aperture radius 
of 20" for the source extraction, since a larger aperture was af- 
fected by the (mainly hard X-ray) background flare, resulting 
in a hardening of the spectrum and a lowering of the signal- 

9 Portable Interactive Multi-Mission Simulator at 
http://heasarc.gsfc.nasa.gov/ Tools/w3pimms.html 
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TABLE 3 

Fluxes, Luminosities, and for the Sample of z > 4 AGNs 



AGN 


AB x 

1450(1+;) A 


M B 


/ v [2500(l+z) A]" 
(mJy) 


log vL v (2500 A) 
(ergs s~') 


Optical 
Ref. 


fb.5-2 keV b 


logL 2 -10keV 
(ergs s~') 




PSS 0121+0347 


18.5 


-28.3 


0.222 + 0.071 


46.94 


1 


17.02±f|9 


45.51 


-1 65+ 003 


SDSS 0210-0018 


19.3 


-27.7 


0.106 + 0.011 


46.71 


2 


14.37^ 


45.27 


'■ J -0.02 


SDSS 0231-0728 


19.3 


-27.9 


0.062 + 0.030 


46.56 


3 


4 - 68 -L29 


45.21 


1 - DZ -0.06 


PSS 0926+3055 


16.7 


-30.1 


1.167 + 0.042 


47.67 


1 


39.04^ 


45.89 


-1 76+ 003 

•■•'"-aoi 


PSS 1326+0743 


17.2 


-29.6 


0.736 + 0.051 


47.46 


1 


27.89+^ 


45.73 


-1 76+ 003 
L - 10 -0.02 


Q 0000-263 


17.5 


-29.3 


0.579 + 0.058 


47.35 


4 


12 59+ ' 75 


45.72 


- 1 70+° 03 


BR 0351-1034 


18.7 C 


-28.2 


0.185 + 0.019 


46.89 


5 


1 1 69+ 2 ' 22 


45.39 


-1 69+ 004 
1 - oy -0.04 


SDSS 1030+0524 


19.7 


-27.8 


0.076 + 0.008 


46.75 


6 


3 nq+ ' 53 


45.18 




BR 2237-0607 


18.3 C 


-28.5 


0.267 + 0.027 


47.08 


5 


8 - 50 -L27 


45.32 


L - -0.03 


SDSS 1306+0356 


19.6 


-27.9 


0.085 + 0.008 


46.76 


6 




45.07 


, 7 ^+0.05 
i - ' J -0.05 



REFERENCES. — (1) Vignali et al. (2003a); (2) Vignali et al. (2001); (3) Vignali et al. (2003b); (4) Schneider, Schmidt, & Gunn (1989); (5) 
Storrie-Lombardi et al. (1996); (6) Brandt et al. (2002) 

"Assuming an ultraviolet (UV) continuum of the form: f v « v~ 79 (Fan et al. 2001). 

b Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10~ 15 ergs cm~ 2 s . 

c Estimated using an empirical relation between the APM R and AB i45o( i+;)A magnitudes and redshift; see § 2 of Kaspi et al. (2000). 



to-noise ratio. For a nominal background level, increasing the 
aperture radius from 20" to 30" should not have resulted in 
any noticeable change in either photon index or normaliza- 
tion (Kirsch et al. 2004). Background regions were taken to 
be at least as large as the source regions; these were annuli, 
when there were no neighboring sources, and circles other- 
wise. The redistribution matrix files (RMFs; which include 
information on the detector gain and energy resolution) and 
the ancillary response files (ARFs; which include informa- 
tion on the effective area of the instrument, filter transmis- 
sion, detector window transmission and efficiency, and any 
additional energy dependent effects) for the spectra were cre- 
ated with the SAS tasks RMFGEN and ARFGEN, respectively. 10 
We note that the two reflection grating spectrometer detectors 
on XMM-Newton do not have sufficient counts to perform a 
high-resolution spectral analysis for any of our sources. 

The spectra were grouped with a minimum of 10-25 pho- 
tons per bin using the task GRPPHA. For each AGN we used 
XSPEC Vl 1.3.0 (Arnaud 1996) to fit the data with the fol- 
lowing models: (i) a power-law model and a Galactic absorp- 
tion component (Dickey & Lockman 1990), which was kept 
fixed during the fit, and (ii) a model similar to the first with an 
added intrinsic (redshifted) neutral absorption component. In 
all the fits we assumed solar abundances (Anders & Grevesse 
1989) and used the PHABS absorption model in XSPEC with 
the Balucinska-Church & McCammon (1992) cross-sections. 
We fitted the data from each EPIC detector alone, jointly fit- 
ted the data of the two MOS detectors, and jointly fitted all 
three EPIC detectors with these models. In general, the five 

10 For more details see the NASA/GSFC calibration memo 
CAL/GEN/92-002 at http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/docs/ 
memos/cal_gen_92_002/cal_gen_92_002.html 



different fitting results were consistent with each other within 
the uncertainties, and the joint fit of all three EPIC detectors 
had the smallest uncertainty on the spectral parameters. The 
best-fit results for the joint fit of the three EPIC detectors are 
summarized in Table[2](first five sources). For all five sources 
the best-fit model is a Galactic-absorbed power law; an F -test 
carried out on the data does not warrant the addition of an in- 
trinsic absorber in any of the spectra. The 90% confidence 
upper limits on intrinsic 2Vh appear in Table |2] The spectra 
of the five AGNs and their best-fit models appear in Fig. \l\ 
Also shown for each AGN is a confidence-contour plot of the 
T-Nu parameter space. No strong systematic residuals are 
seen for any of the fits. Optical fluxes and luminosities, and 
X-ray fluxes and luminosities of the sources computed using 
the best-fit spectral parameters of Table |2 appear in Table [3] 
(first five sources). 



2.3. Serendipitous Sources 

Since the luminous AGNs under study are probably asso- 
ciated with massive potential wells in the earliest collapsed 
large-scale structures, we searched for candidate z > 4 AGNs 
in the XMM-Newton images of three of our sources with 
SDSS imaging (SDSS 0210-0018, SDSS 0231-0728, and 
PSS 1326+0743) to search for AGN clustering on scales of 
«0.1-5Mpc at z ^4-5 (e.g., see Djorgovski et al. 2003; 
Vignali et al. 2003b; Hennawi et al. 2005). We matched 
X-ray and SDSS sources in those fields and utilized the 
Richards et al. (2002) color criteria as well as a pointlike 
morphology criterion to identify «3 X-ray emitting z Si 4 
AGN candidates for future optical spectroscopy. 
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FIG. 2. — 68, 90, and 99% confidence regions for the photon index vs. 
intrinsic column density derived from joint spectral fitting of our sample of 
eight radio-quiet AGNs. Solid (dashed) contours refer to the entire (common) 
energy range of the AGNs (see §@for more details). 



3. ARCHIVAL SAMPLE OF Z > 4 RADIO-QUIET AGNS 

To improve the statistics on X-ray spectral properties, 
we expanded the sample by adding moderate-to-high qual- 
ity data of five radio-quiet z > 4 AGNs, obtained from the 
XMM-Newton and Chandra archives. A minimum of ~200 
(^100) photons per source from XMM-Newton (Chandra) 
observations was considered useful to our analysis. The ob- 
servation log of the four sources observed by XMM-Newton 
appears in Tabled which also gives the references to the orig- 
inal papers presenting the results of those observations. We 
processed the observation data files of these sources follow- 
ing the analysis procedures described in §|2] As in the case of 
our original sample, we were only able to place upper limits 
upon the intrinsic column densities in these sources. The best- 
fit model parameters for the X-ray spectra of these sources are 
listed in Table |2] 

The fifth observation, a 118 ks Chandra/ ACIS-S ob- 
servation of SDSS J130608.26+035626.3 (hereafter 
SDSS 1306+0356) at z=5.99 in 2003 November 29-30, 
was re-analyzed following the same source and background 
extraction parameters as those used in the original paper 
reporting that observation (Schwartz & Virani 2004) applying 
standard CIAO 11 V3.2 routines to the data. The events were 
grouped into bins that had a minimum of 15 counts per bin, 
and the model fit to the data included a constant Galactic 
absorption component (7Vh=2.07x 10 20 cm -2 ), a power-law 
component, and a neutral Fe Ka line at rest-frame 6.4 keV 
with a width of <T=0.1 keV, since there is a feature in the 
spectrum at that location hinting at the presence of such 
a line (see Schwartz & Virani 2004). The best-fit photon 
index is r=1.82^Q 2 j, the flux normalization at 1 keV is 
(1.23±0.24)xl0~ 6 keV cm" 2 s" 1 keV" 1 , and the Fe Ka 
rest-frame equivalent width is LOiJ g keV. 

The optical fluxes and luminosities, and X-ray fluxes and 
luminosities of the five sources, based upon the best-fit spec- 
tral parameters of Table and of the single Chandra obser- 
vation, are listed in Table [3] Plots of the spectra and best-fit 
models, which are generally consistent with our results, ap- 
pear in the original papers cited in Tableland in Schwartz & 
Virani (2004) for the single Chandra observation. 

" Chandra Interactive Analysis of Observations. See 
http://asc.harvard.edu/ciao/ 



4. JOINT SPECTRAL FITTING OF Z > 4 RADIO-QUIET 

AGNS 

To derive the best possible constraints on the average X-ray 
spectral properties of radio-quiet z > 4 AGNs, we jointly fit- 
ted the spectra of the expanded sample of eight radio-quiet 
AGNs (e.g., see § 6 of Vignali et al. 2005). We com- 
bined the spectra from all the detectors of the observations 
of the eight AGNs, 22 datasets in total, and used XSPEC to 
fit their spectra jointly. We fit two different models: (i) a 
power-law that includes both Galactic and intrinsic absorp- 
tion, and (ii) a power-law with Galactic absorption including 
a Compton-reflection component, using the PEXRAV model in 
XSPEC (Magdziarz & Zdziarski 1995), and a neutral narrow 
Gaussian Fe Ka line. The Compton-reflection component, 
also known as a Compton-reflection 'hump', is the spectral 
manifestation of hard X-ray photons emitted in a corona of 
hot electrons and reflected off the relatively colder accretion 
disk (e.g., see § 3.5 of Reynolds & Nowak 2003 and refer- 
ences therein). The Compton-reflection hump feature is ob- 
served in the rest-frame ^7-60 keV energy range, peaking 
at rest-frame ^30 keV. The joint-fitting process was carried 
out twice for both models; first, to include the entire energy 
range of all the datasets, and a second time to include only 
the common rest-frame energy range of the eight sources: 
1.5 ^ iirest ~ 51 keV. The total number of photons used in the 
joint-fitting process is ^7000 (^6400 for the common energy 
range), and it is larger by an order of magnitude than the num- 
ber of counts used in the joint fitting of 48 z > 4 radio-quiet 
AGNs by Vignali et al. (2005). 

The fitting process with the first model required that all 
datasets be fitted with a single power-law and a single in- 
trinsic absorption column. Therefore the photon index and 
intrinsic absorption column were each tied to a single value 
across all data sets, and these two parameters were allowe d to 
vary jointly (see the discussion of intrinsic dispersion in $ 15.11 
as partial justification of this approach). The power-law nor- 
malizations of all data sets were allowed to vary freely. For 
the second model we used a single power-law and a single 
relative -reflection component (R, where R is defined as Q./27Z 
and £2 is the solid angle subtended by the reflecting medium) 
of the Compton hump for all datasets, and allowed them to 
vary jointly. We fixed the rest-frame energy and width (c) 
of the neutral Fe Ka line to 6.4 and 0. 1 keV, respectively, and 
tied all the power-law and Gaussian normalizations to two dif- 
ferent values, allowing them to vary jointly; each dataset was 
assigned a constant scaling factor to account for the different 
flux levels of each AGN and each detector. The best-fit val- 
ues and upper limits for the parameters used in the joint fits 
are given in Table|4] Fig.|2]is a contour plot of F versus Nh 
as a result of the joint-fitting process using the first model. 
The noticeable difference between the confidence contours in 
the entire and common energy ranges stems from the fact that 
the two AGNs at z —6 increase the low-energy threshold in 
the common energy range, and as a consequence many low- 
energy photons from the X-ray bright AGNs at z ~4.2 are 
ignored in the fit. This causes the upper limit on the intrinsic 
absorption in the common energy range to rise significantly 
relative to the upper limit when using the entire energy range. 

We repeated the joint-fitting process with the first model 
and fitted the four most optically luminous radio-quiet sources 
separately from the four least optically luminous radio-quiet 
sources both for the entire and common energy ranges of 
the two sub-samples. Both F and Nn in the two luminosity 
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TABLE 4 

Best-Fit Parameters from Joint Fitting of Eight Radio-Quiet z > 4 AGNs 



Spectral Model 
Galactic-Absorbed Power Law and 


Energy 
Range 3 


N H 

(10 22 cm" 2 ) 


r 


EW(Fe Ka) b 
(eV) 




Z 2 (DOF) 


Intrinsic Absorption 


E 


<0.28 


, Q7 +0.06 






308.3 (360) 


Intrinsic Absorption 


C 


<0.84 


2 02+ 04 

z ' uz -0.06 






291.3 (338) 


Compton-Reflection and neutral Fe Ka line 


E 




1 98 +0J0 

1,y °-0.09 


<189 


< 1.22 


306.0 (359) 


Compton-Reflection and neutral Fe Ka line 


C 




1 nc+0.10 


< 191 


< 1.03 


291.2 (337) 



a E - entire energy range of all the spectra. C - the common rest-frame energy range among all sources (1.5-51 keV). 

b Rest-frame equivalent width calculated with a mean redshift of 4.8. 

c Compton hump relative reflection parameter (see Magdziarz & Zdziarski 1995). 



groups are consistent with one another and with the values 
we obtain by fitting the entire radio-quiet sample. For the fit- 
ting using the entire energy range, we measure r=1.99±0.06 
for the most-luminous sources and r=1.95^Q2o for the least- 
luminous sources. 

5. RESULTS AND DISCUSSION 
5.1. X-ray Photon Index 

The photon indices and their uncertainties for the full sam- 
ple of ten AGNs, in the ^1-50 keV rest-frame energy range 
(the range where most photons are detected), are listed in 
Table [2] A statistical summary of our expanded sample's 
photon-index distribution appears in Table [5] Using the max- 
imum likelihood method of Maccacaro et al. (1988) we find 
no significant intrinsic dispersion in T; this is expected since 
the individual uncertainties on Y are larger than the spread 
in individual Y values. In Fig. [3] we plot Y versus abso- 
lute B magnitude and redshift for our expanded sample, to- 
gether with previous measurements of the photon index in 
other AGN samples (see § 6 of Vignali et al. 2005 and ref- 
erences therein). There is no obvious systematic difference 
between the photon indices of our radio-quiet sources and 
those of luminous AGNs observed at lower redshifts (e.g., 
the Piconcelli et al. 2005 sample of nearby AGNs). Using 
a Spearman rank correlation method we find no significant 
correlations between Y in radio-quiet sources and either lu- 
minosity or redshift as also found at lower redshift (z ^ 0.4) 
by, e.g., Porquet et al. (2004). Applying a linear regression 
algorithm that accounts for intrinsic scatter and measurement 
errors (Akritas & Bershady 1996), we obtain \dY/dz\ < 0.04 
for the 51 radio-quiet AGNs plotted in Fig. [3] From the joint 
fitting of eight radio-quiet AGN spectra using the absorbed 
power-law model we find Y = 1.97+QQ4. This is consistent 
with both the unweighted mean of the individual Y values and 
the counts-weighted mean Y for the radio-quiet sample. 

The mean photon indices we find for z > 4 radio-quiet 
AGNs are also all consistent with those calculated by Vig- 
nali et al. (2005). All these photon-index values confirm and 
strengthen the Vignali et al. (2005) result that the X-ray pho- 
ton index does not depend on optical luminosity or redshift. In 
particular, our results do not support either of the reports for 
a significantly different photon index at low versus high red- 
shift (Bechtold et al. 2003; Grupe et al. 2005). Bechtold et al. 
(2003) used a heterogeneous sample of low-redshift ROSAT 
sources and high-redshift sources detected with Chandra. The 



mix of different datasets may have produced a bias which led 
Bechtold et al. (2003) to report a flattening of Y at higher 
redshift (e.g., due to spectral curvature at low energies and 
mission-to-mission cross-calibration uncertainties). In ad- 
dition, the counts-weighted mean Y of their z > 4 Chandra 
sources is 1.50±0.15 (see their Table 2). This is signifi- 
cantly lower than the average T of 1 .93 in 09 f° r ^8 radio-quiet 
AGNs, including the Bechtold et al. (2003) sources, obtained 
by Vignali et al. (2005). Comparing the photon indices of our 
four sources which overlap with those of Grupe et al. (2005) 
shows that in three cases the values are in good agreement. 
However, for SDSS 1030+0524 the T=2.65±0.34 found by 
Grupe et al. (2005) is significantly steeper, and is inconsis- 
tent both with the value we obtain, r=1.92+02j, and that ob- 
tained by Farrah et al. (2004), r=2. 12+0. 11. The photon in- 
dex we obtain for SDSS 1030+0524, while slightly harder, 
is compatible within the errors with the one found by Far- 
rah et al. (2004). By inspection of Tables 1 and 4 of Grupe 
et al. (2005) it is also apparent that the relatively steep mean 
r they report for their radio-quiet sources, 2.23+0.48, is an 
unweighted mean dominated by sources with limited photon 
statistics («50 counts), which in most cases result in best-fit 
r exceeding a value of 2.2. 

The insignificant dispersion in individual Y values for our 
z > 4 AGN sample may partially be driven by a selection bias, 
since we specifically targeted some of the most X-ray lumi- 
nous AGNs known. However, observing down the AGN lu- 
minosity function at z > 4 is challenging with existing X-ray 
observatories. In fact, the faintest X-ray sources at those red- 
shifts with enough counts to allow meaningful measurement 
of the photon index were detected in the 2 Ms Chandra Deep 
Field North survey (Vignali et al. 2002; see Table 1 of Brandt 
et al. 2005 for a list of X-ray faint AGNs at z > 4). The best 
studied source in that survey, CXOHDFN J123647.9+620941 
at z=5. 19 and with L 2 -io keV=10 44 ergs s -1 , has r=1.8±0.3, 
which is consistent with the photon indices we find for our 
sources (Vignali et al. 2002). The photon indices of our least 
luminous sources do not appear to be significantly different 
from those of the more luminous ones both individually and 
on average by splitting the joint spectral fit into two lumi- 
nosity groups (see §0}. Inspection of Table 4 of Vignali et al. 
(2005) also does not show a strong dependence of Y on optical 
luminosity. In spite of being extremely luminous, our z > 4 
AGNs are otherwise representative of the AGN population by 
having 'normal' rest-frame UV properties (i.e., no broad ab- 
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FIG. 3. — The X-ray photon index (T), versus (a) absolute B magnitude and (b) redshift, adapted from Vignali et al. (2005). The eight (two) radio-quiet 
(moderately radio-loud) AGNs presented in this work are marked with open (filled) squares. Note that the photon index does not show any clear dependence on 
either luminosity or redshift. 



TABLE 5 

Photon Index Statistics of z > 4 AGNs 



Photon Index (Y) 



Statistic 


Entire Sample 


Radio-Quiet Sources 


Unweighted Mean 


1.89 


1.92 


Unweighted Standard Deviation 


0.11 


0.11 


Unweighted Error of the Mean 


0.03 


0.04 


Counts-Weighted Mean 


1.93 


1.95 


Median 


1.86 


1.90 



sorption lines, typical emission-line properties and rest-frame 
UV continua; Vignali et al. 2001, 2003a, b). Recent multi- 
wavelength observations of z > 4 AGNs have found that their 
overall spectral energy distributions (SEDs) are not signifi- 
cantly different from those of lower redshift sources implying 
no SED evolution (e.g., see Carilli et al. 2001 and Petric et al. 
2003 for radio observations; Pentericci et al. 2003 and Van- 
den Berk et al. 2001 for UV-optical observations). We argue 
that the broad-band X-ray spectral properties of z > 4 AGNs 
are also not significantly different than those of lower redshift 
sources. 

5.2. Intrinsic Absorption 

In all the individual fits of the z > 4 AGNs we obtained 
an intrinsic Nh consistent with zero with 90% confidence. 
Table |2] lists the upper limits on the intrinsic absorption; 
they are all in the range «10 22 -10 23 cm~ 2 . The joint-fitting 
procedure yielded an upper limit on the common intrinsic 
Nh of ~(3-8)x10 21 cirr 2 (Table |U and Fig. EJ. These con- 
straints on intrinsic absorption are the tightest presented to 
date for z > 4 AGNs, and our individual-fitting results are 
less prone to bias than the joint-fitting results that have been 
predominantly utilized previously. 12 Our results show that 
radio-quiet AGNs at z > 4 do not appear to have higher ab- 
sorption columns than their lower redshift counterparts. The 
metallicity in the inner regions of high-redshift AGNs may be 
higher, on average, than in their lower-redshift counterparts 

12 Since photoelectric absorption by neutral gas is an extremely strong 
spectral effect, a column density derived via joint fitting of several sources 
may not always be a good estimate of their 'typical' column density. The 
presence of even one source with significant low-energy flux can bias the 
joint-fitting column density downward, since such flux is highly inconsistent 
with a strong photoelectric-absorption cutoff at low energy. 



mainly due to their higher luminosities (e.g., Dietrich et al. 
2003). Since X-ray photoelectric absorption is predominantly 
due to metals and not hydrogen, the inferred intrinsic neutral 
hydrogen column density depends upon the assumed metal- 
licity. Therefore, Nh values are inversely proportional to the 
assumed metallicity, and the quoted upper limits on the best- 
fit intrinsic columns given in Tableware conservative values 
(as solar abundances were assumed). 

Despite the considerable large-scale evolution the Universe 
has experienced since z ~6, and in particular the strong evo- 
lution in the AGN population (e.g., Fan et al. 2003; Croom 
et al. 2004; Barger et al. 2005; Brandt & Hasinger 2005), we 
do not detect significant changes in the central regions of lu- 
minous radio-quiet AGNs. Since the host galaxies of AGNs 
at high redshifts are probably still in the formation process, 
one might expect large amounts of gas and dust in the inter- 
stellar medium of these hosts that could increase the amount 
of line-of-sight absorption. Signatures of large amounts of in- 
terstellar molecular gas and dust have been recently detected 
in some of the most distant AGNs (e.g., Carilli et al. 2001; 
Bertoldi et al. 2003; Walter et al. 2003; Maiolino et al. 2004). 
In addition, the fact that radio-loud AGNs show signs for in- 
creased intrinsic absorption with redshift (see § Q), makes it 
of interest to check if there is a similar trend for the radio- 
quiet population. Our results do not show significant intrinsic 
absorption in any of the z > 4 radio-quiet AGNs, which is in 
broad agreement with, e.g., the Laor et al. (1997) results for 
local radio-quiet AGNs. We caution that our AGN sample is 
far from being complete or representative of the entire AGN 
population. We have specifically considered only optically 
and X-ray bright, radio-quiet AGNs with no broad absorption 
lines. Steffen et al. (2003) have found that the fraction of type 
1 (i.e., unobscured) AGNs increases with X-ray luminosity, 
suggesting luminosity-dependent absorption. It is therefore 
possible that the relatively lower intrinsic absorption in our 
sources is a consequence of them being among the most lu- 
minous AGNs known. 

5.3. Fe Ka Line and Compton Reflection 

The joint-fitting process placed respectable constraints on 
the strength of a putative neutral narrow Fe Ka line; a 
90% confidence upper limit on the rest-frame equivalent 
width is ^190 eV. We also constrained the relative reflec- 
tion component of a Compton-reflection hump (see Table|4j. 
These are the best constraints to date on these properties 
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o Stratcva et al. (2005) 
□ This work (radio quiet) 
■ This work (radio loud) 



28 29 30 31 32 

Log L v (2500A) [ergs s'' Hz''] 

FIG. 4. — The optical-to-X-ray SED parameter, o^x, versus luminosity 
density at 2500 A. The ten AGNs studied in this work (squares) follow the 
Strateva et al. (2005) a ox -L v (2500 A) relation (solid line), a m = -0.136 x 
Lv(2500 A)+2.616. The two moderately radio-loud sources presented in this 
work, PSS 0121+0347 and SDSS 0210-0018, are the only such sources on 
this diagram and are marked with filled squares. The reason why most of 
our sources lie above the best-fit line may be attributed to the fact that we 
specifically targeted X-ray bright sources. 

of z > 4 AGNs. Due to the high mean X-ray luminosity 
((Log L2-10 kev(ergs s -1 )) ~45.5) of our sample, both the 
Fe Ka line and Compton hump are expected to be relatively 
weak. Using the recent Page et al. (2004a) anticorrelation 
between EW(Fe Ka) and i2-i0kev (the so-called Fe Ka 
"Baldwin Relation"), we find an empirically expected range 
of EW(Fe Ka)~25-80 eV for our sources, several times 
lower than our upper limit (but see also Jimenez-Bailon et 
al. 2005). Given our results, the EW(Fe Ka)«l keV for 
SDSS 1306+0356 (see §|5J appears remarkable and should 
be investigated with better data. The constraints we ob- 
tain on the Compton-refiection component, R 1.2, are con- 
sistent with those obtained for some z ~0-l AGNs (e.g., 
Page et al. 2004b), but are significantly lower than the value 
fl=2.87±0.96 found in PG 1247+267 at z=2.04, the most lu- 
minous source in which such a component was significantly 
detected (Page et al. 2004b). 

5.4. Optical-to-X-ray SED 

For each of the ten AGNs we computed the optical-to-X-ray 
SED parameter, a ox , defined as: 

„ _ lo g(/2kev// 2 500 A) m 
Lt x — if / \ {J 

log(V2kev/v 25 o oA ) 

where fa k e v and /2500 a are me f ux densities at rest-frame 
2 keV and 2500 A, respectively (e.g., Tananbaum et al. 1979). 

The Oox values and their Iff errors appear in Table [5] 
The errors on Oox were computed following the "numerical 
method" described in § 1.7.3 of Lyons (1991), taking into 
account the uncertainties on the X-ray power-law normaliza- 
tions and photon indices, and the effects of possible changes 
in the UV-optical slope (from a=-0.5 to -0.79). In Fig. |4] 
we plot Oox for our sources against their optical luminosities, 
and include a sample of 228 radio-quiet AGNs (mainly from 
the SDSS) from Strateva et al. (2005) in the diagram. Also 
plotted in Fig.0|is the best-fit line of the Strateva et al. (2005) 
a ox -M2500 A) relation (Eq. 6 of Strateva et al. 2005). The 
a ox values of the ten AGNs presented in this work, which 
represent extreme luminosities and redshifts, generally follow 
that relation, although most of them lie above the best-fit line. 
This may be due to a selection bias since we have specifically 
targeted X-ray bright sources. 




Log F ft 5 2 keV 1 er 8 s cm 2 s ! ] First Epoch 



FIG. 5. — Two-epoch Galactic-absorption corrected 0.5-2 keV fluxes for 
seven of the z > 4 AGNs in our sample. The solid line marks the 1:1 flux 
ratio, and the two dotted lines mark 1:2, and 2:1 flux ratios, to guide the eye. 
SDSS 0231—0728 clearly varied by more than a factor of two between the 
two epochs. The second and third most variable sources, PSS 0121+0347 
and SDSS 0210—0018, are moderately radio loud and are marked with filled 
circles. 

5.5. X-ray Variability 

Using the event files of our XMM-Newton observations, 
we searched for rapid (~ 1 hr timescale in the rest frame) vari- 
ability of the AGNs by applying Kolmogorov-Smirnov tests 
to the lists of photon arrival times. No significant variations 
were detected, and the fractional variability of the count-rate 
light curves was smaller than 15% in all cases. To test if the 
X-ray flux exhibits long-term (months-to-years) variations in 
our sample, we list in Table [6]pairs of two-epoch 0.5-2 keV 
fluxes for seven z > 4 sources which were observed with 
Chandra in the first epoch, and with either XMM-Newton or 
Chandra in the second epoch. The relative uncertainties be- 
tween flux measurements obtained with XMM-Newton ver- 
sus those obtained with Chandra are considered to be the 
smallest among different X-ray observatories, and are nor- 
mally < 10% (Snowden 2002; S. Snowden, 2005, private 
communication). Table [6] also lists the rest-frame temporal 
separation between the two observations, and a % 2 statistic to 
quantify the significance of variation between the two epochs. 
The null hypothesis is that the flux in each epoch is equal to 
the mean flux of the two epochs. A larger % implies higher 
significance of variability, and values larger than 2.71 (6.63) 
imply that the two fluxes are inconsistent within their uncer- 
tainties, i.e., there is significant variability at a 90% (99%) 
confidence level (e.g., Avni 1976). Of the seven AGNs studied 
only two sources, PSS 1326+0743 and SDSS 1030+0524, do 
not show a significant (i.e., at >90% confidence) long-term 
variation. Fig.[5]plots the fluxes of the second epoch against 
those of the first, as given in Table |6] 

While most AGNs varied by no more than a factor of «2 
between the two epochs one source, SDSS 023 1 —0728, faded 
by a factor of ^4 between the first {Chandra) observation and 
the second (XMM-Newton) one. This flux change occurred 
over a rest-frame period of 73 d. This is the largest change 
in X-ray flux observed in a z > 4 radio-quiet AGN (e.g., Pao- 
lillo et al. 2004). Given the optical flux of the source, and 
using the Vignali et al. (2003b) relation between optical and 
X-ray fluxes, it is Likely that this source was caught in an 
X-ray high state in the first epoch (Vignali et al. 2003b), since 
its X-ray flux in the second epoch (this work) agrees with 
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TABLE 6 

Two-Epoch X-ray Fluxes Obtained by XMM-Newton and Chandra 



AGN 


First Epoch 

Fo.5-2 keV a 
(10~ 15 ergs cm~ 2 s _1 ) 


Ref. 


Second Epoch 

Fo.5-2 keV a 
(10~ 15 ergs cm -2 s - ') 


Ref. 


Af b 
(days) 


x 2c 


PSS 0121+0347 


34.7jft 


1 


17.02l 2 ^ 


2 


135 


20.40 


SDSS 0210-0018 


28.11JJ 


3 


14 37 +L77 


2 


258 


16.35 


SDSS 0231-0728 


zl -'-4.6 


4 


4.6811" 


2 


73 


41.42 


PSS 0926+3055 


27.71J] 


1 


39.04±|g 


2 


186 


6.43 


PSS 1326+0743 


24 4+ 4 ' 3 


1 


27.891';* 


2 


139 


0.64 


SDSS 1030+0524 




5 


3 09 +a53 


6 


68 


0.45 


SDSS 1306+0356 


4 63+177 


5 


2 7 4+0.56 


7 


96 


3.35 



REFERENCES. — (1) Vignali et al. (2003a); (2) This work; (3) Vignali et al. (2001); (4) Vignali et al. (2003b); (5) Brandt et al. (2002); (6) 
Farrah et al. (2004); (7) Schwartz & Virani (2004) 

"Galactic absorption-corrected flux in the observed 0.5-2 keV band. Errors represent 90% confidence limits on the fluxes. 
b Rest-frame time differen ce be tween the two epochs. 
Variability statistic; see § 15.51 



the value predicted from its optical flux (assuming the optical 
flux is nearly constant). Vignali et al. (2003b) also noted that 
SDSS 023 1 —0728 was X-ray brighter than expected (see their 
Fig. 5). The spectral slope of the source also shows a possible 
indication of flattening from r=2.8^ ^ to r=1.85+g ^ be- 
tween the two epochs, but the significance is only ^1(7 due to 
the limited number of counts (^25) in the first Chandra snap- 
shot observation. This is a tentative indication for a transition 
from a soft/high state to a hard/low state in this source, as has 
been seen for a few local AGNs (e.g., Guainazzi et al. 1998; 
Maccarone, Gallo, & Fender 2003). 

6. SUMMARY AND FUTURE PROSPECTS 

Moderate-to-high quality X-ray spectra of ten luminous 
AGNs at z > 4, including five new XMM-Newton observa- 
tions, have been analyzed to extract X-ray spectral proper- 
ties and search for flux variations. We found that the X-ray 
power-law photon index does not depend on either luminos- 
ity or redshift, broadly consistent with the lack of multiwave- 
length AGN SED evolution, and it has no detectable intrin- 
sic dispersion at z > 4. Upper limits placed on the intrinsic 
absorption in eight radio-quiet z > 4 AGNs show that, on av- 
erage, the X-ray light of luminous radio-quiet AGNs at z > 4 
is not more absorbed than in nearby AGNs. All this suggests 
that the X-ray production mechanism and the central environ- 
ment in radio-quiet AGNs do not depend on luminosity and 
have not evolved over cosmic time. We also place constraints 
on the strength of a putative narrow Fe Ka line, and constrain 
the relative strength of the Compton-reflection hump at high 
rest-frame X-ray energies. Tighter constraints on those X-ray 
properties of luminous AGNs at z > 4, with existing X-ray ob- 
servatories, can only be achieved by lengthy (several 100 ks) 
observations of the most luminous sources. 

A search for short-timescale (~1 hr) X-ray variability found 
no significant variations, but on longer timescales (months- 
to-years) we have found significant variations in five of our 
sources, noting one extreme radio-quiet source in which the 
X-ray flux dropped by a factor of ~4 over a 73 d rest-frame 



period; this is the most extreme X-ray variation ever ob- 
served for a radio-quiet AGN at z > 4. X-ray monitoring of 
SDSS 0231-0728 and other AGNs at high redshift will en- 
able us to link the variability properties of high-redshift AGNs 
with some of their fundamental properties, such as BH mass, 
luminosity, and accretion rate, as is done at lower redshift 
(e.g., O'Neill et al. 2005). Two recent studies have reported 
that the X-ray flux variations of high-redshift AGNs are larger 
than those of their nearby counterparts (Manners et al. 2002; 
Paolillo et al. 2004). The data presented in this paper do not 
allow us to test this claim. Multiple-epoch X-ray observations 
of a large sample of high-redshift radio-quiet AGNs are re- 
quired to investigate their variability properties, and compare 
them with those of the local AGN population. 
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